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Abstract. We study domino tilings on the Aztec diamond, with two different weights for 
horizontal and vertical dominos. Using a result of Helfgott, and a one-to-two lifting from 
the correlation kernel of a related interlaced particle system, we guess and then prove a 
formula for the inverse of the Kasteleyn matrix of the Aztec diamond. This lets us compute 
asymptotics for the determinantal southern domino process, using Kenyon's local statistics 
formula. We find that at the northern edge of the chaotic region, the southern domino 
process converge to a thinned version of the Airy point process. Conversely, at the southern 
edge of the chaotic region, the holes of the southern domino process converge to a thickened 
Airy point process, where one adds multiple points to the Airy point process. We also 
outline the proof that the fiuctuations of the height function converge to the Gaussian free 
field. 



1. Introduction 

The Aztec Diamond of order n is a planar region which can be completely tiled with 
dominos — two-by-one rectangles. Over the past twenty years, this particular shape has 
come to occupy a central place in the literature of domino tilings of plane regions. Tilings 
of large Aztec diamonds exhibit striking features — the main one being that these tilings 
exhibit a limit shape, described by the so-called Arctic circle theorem |17j . The object of 
study in this paper is a perturbation of the uniform measure on domino tilings, in which 
vertical tiles are preferred over horizontal tiles. See Figure [l] for pictures of tilings of a 
relatively large Aztec diamond. 

There are several alternate descriptions of a tiling of an Aztec diamond. A domino tiling 
is equivalent to a perfect matching, or dimer cover, on the dual graph of the region which 
is tiled. There is an equivalent family of nonintersecting lattice paths, called DR paths [32], 
and there is a description as a stack of a certain sort of blocks, called Levitov blocks \30\ [31] 
(see Figure [T]) There is also a well-studied interlacing particle process which is equivalent to 
the tiling model in a certain sense, but this equivalence is not bijective: there are 2"(""'"^)/^ 
tilings of the Aztec diamond, whereas the number of configurations of the particle process 
are equinumerous with order-n Alternating sign matrices or configurations of the six-vertex 
model. However, the correspondence is weight preserving and locally defined; it maps certain 
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collection of tilings to a configuration of the free-fermion six- vertex model [13] , preserving 
the relative weights. 

It is this point process whose local asymptotics have been studied most thoroughly [19]; 
in particular, the boundary of the frozen region is more easily described using the particles, 
since it is related to the position of the last particles on different lines. There is also a 
relationship between these particles and a certain sort of zig-zag path in the tiling (distinct 
from the zig-zag paths studied in [7]) which is helpful. However, the many-to-one nature of 
the correspondence necessarily loses some of the information about the original tilings. In 
both the domino and particle pictures, we get determinantal point processes, although the 
precise combinatorial relation between the two processes is not immediate. 

We correct this situation in this paper, by giving a formula for the inverse Kasteleyn 
matrix for the Aztec diamond. This generalizes a previous result by Helfgott [16] to the case 
when the horizontal and vertical dominos have different weights. Using an observation of 
Kenyon [26] , it is then possible to compute probabilities for various configurations of dominos 
and their asymptotic limits as the size of the Aztec diamond increases. In particular, we 
answer two questions about domino tilings of the Aztec diamond which would be much 
harder to state, or prove, in the interlaced particle process: 

(1) How does the boundary of the frozen region manifest itself if we can only see one 
type of dominos, for example, the southern dominos? 

(2) Do the fluctuations of the chaotic region of the Aztec diamond tend, as one would 
expect for any dimer model, to the Gaussian Free field? 

1.1. The southern domino process: North Boundary. There are in fact four different 
types of dominos in a tiling: the dominos can be placed in two orientations, each of which 
comes in two different parities (determined by the bipartition of the dual graph on which 
the dominos are placed). Due to the Arctic circle theorem, with probability one, there are 
only dominos of one of these four types clustered near each of the four corners of the Aztec 
diamond. For this reason (and others, see [13]) we call the four types of dominos north, 
south, east, and west, which are coloured red, green, yellow and blue in Figure [T] in the 
electronic version of this article. 

For the moment, ignore all but the south dominos (the green ones in Figure [T]). Viewing 
each domino as a point, the set of all southern dominos form a determinantal point process. 
Note that the positions of the southern dominos do not specify the tiling uniquely so they 
only give a partial description of the tiling (though, together with any of the other types of 
dominos, they do). We analyze the distribution of the southern dominos along a diagonal 
line in a large Aztec diamond, scaled so as to study the two intersections between this line 
and the frozen boundary of the tiling (an "arctic ellipse" , in the weighted case. It would 
be possible to extend what we have done to analyze the joint distribution of all southern 
dominos in the Aztec diamond; as is, our analysis extends previous results in [8] on placement 
probabilities of single dominos. 

We will show that the appropriate scaling limit of the point process of southern dominos 
along a diagonal line close to the boundary of the northern frozen region is given by a 
thinned Airy kernel point process, where the amount of thinning depends on the relative 
weight of the horizontal and vertical dominos. This can be heuristically understood in the 
following way. The Airy kernel point process, as mentioned earlier, is the edge limit of the 
particle process along a diagonal line, and these in turn are given by the intersections of the 
non- intersecting paths and the diagonal line. Sometimes these intersections occur along a 
southern domino, and sometimes not. If we only see the southern dominos we only see some 
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of these intersections and which of them we see is essentially random; so we might expect, 
in the limit, a random thinning of the Airy kernel point process. A priori it is not clear that 
the thinning becomes independent in the limit, but this turns out to be the case. 

1.2. The southern domino process: South Boundary. If, instead, we examine the 
southern frozen region, we find that almost all of the dominos are coming from the south 
domino process, and thus lie in a frozen, brickwork pattern predominantly. The southern 
boundary is a "hole" in this regular pattern. Consider the holes between the south dominos 
along a diagonal line in a neighborhood of the southern boundary. These holes also form a 
determinantal point process but in a scaling limit it does not converge to a simple process, 
but rather to a multiple point process with independent geometric probabilities for the points 
in an Airy kernel point process. This can be seen in Figure[Tj there is a tendency for dominos 
of like types to cluster together along the southern boundary. This tendency continues even 
in the limit, with a cluster of k dominos becoming a point of multiplicity k. The multiplicity 
increases as we go towards the lower tangency point of the arctic ellipse, a fact which can 
also be observed in Figure [TJ 

We will also consider some scaling limits when the relative weight of vertical and horizontal 
dominos varies with n. Furthermore we will use our knowledge of the inverse Kasteleyn 
matrix to show the convergence of the height function to a Gaussian Free Field. 

1.3. Previous work. Domino tilings on the Aztec diamond were originally introduced 
in [13] as a model connected with the alternating sign matrices. In this section, we give 
only a partial overview of the literature on the asymptotics of domino tilings of the Aztec 
diamond. 

Due to the existence of limit shape of large tilings [HI [19], the Aztec diamond was one of 
the main motivating examples for the general theory of the existence of a limit shape for 
general tiling models on bipartite graphs [HJ [U [29]. For these particular tiling models, one 
can define a height function which is in bijection (up to height level) with the tiling ^38j . The 
behavior in the bulk holds from the general results in |29j which says that the measure of 
domino tilings is locally a Gibbs measure and dependent on the average slope of the height 
function at the specified point. 

The edge behavior, that is, the behavior between the frozen and unfrozen regions has, 
been of particular interest to the random matrix community, as the fluctuations are of size 
n^^^. This is the same size as the fluctuations of the largest eigenvalue of the Gaussian 
Unitary Ensemble (GUE). Indeed, |18| 119] showed that the law of the particles associated 
to the tiling is given by the Airy process and that the position of the last particle is given 
by the Tracy Widom distribution, F2, see for example [Ij. Furthermore, |2l| showed that 
the distribution of these particles becomes the GUE minor process at the intersection of the 
liquid region and the boundary of the Aztec diamond while |i4j obtained the similar results 
for a certain half Aztec diamond. 

There are a handful of other explicitly inverted Kasteleyn matrices in the literature. The 
inverse of the Kasteleyn matrix of the Aztec diamond was computed by Helfgott [T6] in the 
case of the uniform measure on domino tilings; the results in [^ [2^ rely on explicit inverses 
of four Kasteleyn-like matrices that, together, count perfect matchings on a torus-embedded 
graph. Finally, Kasteleyn computes the eigenvalues and eigenvectors of the Kasteleyn 
matrix of the m x n grid graph explicitly. 
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Figure 1. Left a tiling of an Aztec diamond with a = ^, and the corre- 
sponding dimer cover. The green dominos along each horizontal row give the 
southern domino process. Right: the height function associated to this tiling 
(realized pile of Levitov blocks \30\ l3T] ) . 
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2. Preliminaries and Results 

2.1. Determinantal Process. Determinantal point processes are a key part of the analysis 
used in this paper. Here, we briefly describe these processes but more in depth treatises of 
determinantal point processes can be found in |20j and [37j. 

Let A be a Polish space and take A^(A) to denote space of counting measures ^ ovl E 
with ^(-B) < oo for every boundary i? C A. A point process on A is a probability measure 
P on A^(A). Let M„ denote the factorial moment measure, that is, for disjoint Borel sets 
Ax,..., Am in A, 

(2.1) M„(A^i X ••• X A;^") =E 

Suppose that A is a reference measure on A. For example, if A = M we can choose A to 
be the Lebesgue measure. If 

(2.2) Mn{Ai, . . . , An) = pn{xi, . . . ,Xn)dX{xi) . . .dX{Xn) 

J AiX---xA„ 

for all Borel sets Ai in A, we call pn to be the n*'* correlation function. For discrete pro- 
cesses, Pnixi, . . . , Xn) is equal to the probability of seeing particles at . . . x„, whereas for 
continuous processes, pn is the density of seeing particles. For example, if Pn{xi, . . . ,x„) = 
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p(xi) . . . p{xn) where p £ with A = M and A is the Lebesgue measure, then the point 
process is the Poisson point process on M. 

A point process is called determinantal if there exists a function K : A x A — )• C called 
the correlation kernel, with 

(2.3) Pn{xi, . . . , x„) = det (K{xi, Xj))'^ .^^ . 

This leads to the following characterization of a determinantal point process. Let C^(A) be 
the set of all non-negative continuous functions on A with compact support. Take ip G C^{A), 
let A denote the support of ijj and set (j) = 1 — e~^. Let 1a denote the indicator function for 
the set A. Then, provided (j)KlA is trace class and the Fredholm determinant is given by its 
Fredholm expansion, 

(2.4) E [e-^. '^^^^)] = det(/ - cPKlA)mA,x), 
where xj are the points in the process. 

2.2. Kasteleyn Matrix for the Aztec Diamond. For a bipartite graph, the Kasteleyn 
matrix, introduced in [23j to count perfect matchings of a graph, is a signed weighted (possi- 
bly complex entried) adjacency matrix with rows indexed by the black vertices and columns 
indexed by the black vertices. We describe the Kasteleyn matrix for the Aztec diamond 
below but first remind the reader of two theorems. 

Suppose that K denotes the Kasteleyn matrix for a finite bipartite graph G. 

Theorem 2.1 ([23]). The number of weighted dimer coverings of G is equal to \ detK\ 

Suppose that E = {ei}"^| are a collection of distinct edges with = {bi,Wi), where hi 
and Wi denote black and white vertices. 

Theorem 2.2. The dimers form a determinantal point process on the edges of G with 
correlation kernel given by 

(2.5) L{ei,ej) = K{h,Wi)K'\wj,bi) 
where K(b,w) = Kbw and K^^{w,b) = {K~^)wi, 

Proof. By [24J, we have that the probability of finding dimers at the edges ei, . . . , e„ is 

n 

P(ei, . . . , = n ^(^^' ^0 {K-Hwj, bi))l.^^ 

(2, 6l 

= det {K{b„w,)K-\vOj,b,))l^^^ 
= det {L{ei,ej))l.^-^ 

□ 

An Aztec diamond of size n is a region of lattice squares described by {(x, y) : -|- |y| < 
n + 1}. A domino is a two by one rectangle that can be placed vertically or horizontally 
on the Aztec diamond. A domino tiling of the Aztec diamond is a covering of the Aztec 
diamond so that each square is covered exactly once by a domino. An equivalent viewpoint 
is to replace the Aztec diamond by its dual graph and a domino replaced by a dimer which 
is an edge. We rotate the Aztec diamond by 7r/4 counter clockwise. As the Aztec diamond, 
its dual graph, a domino and a dimer are synonymous, we will interchange between the 
definitions. Figure [2] shows the rotated Aztec diamond and a dimer covering of the dual 
graph. 
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Figure 2. The figure on the left shows an Aztec diamond of size 3 with a 
dimer covering of the dual graph. The domino tiling can be seen by placing 
dominos over the dimers. The figure on the right shows the co-ordinates of 
the Aztec diamond (dual graph) - the whtie vertices x = (xi,X2) have xi 
mod 2 = 1 and X2 mod 2 = while the black vertices are the remaining 
vertices. 

We give the (dual graph) Aztec diamond of size n the following co-ordinates: let W and B 
denote the set of white and black vertices for an Aztec diamond of size n. We have that 

(2.7) W = {(xi, X2) : xi e 2Z + 1, X2 G 2Z, 1 < xi < 2n - 1, < X2 < 2n} 
and 

(2.8) B = {(1/1,2/2) : yi G 2Z, ?/2 G 2Z + 1, < yi < 2n, 1 < y2 < 2n - 1}. 

Set ei = (1, 1) and 62 = (—1, 1). We have that x G W is adjacent toyGBiffy — x = ztcj for 
i G {1,2}. See Figure [2] for the co-ordinates of an Aztec diamond of size 3. 

For w G W and 6 G B, a dimer {w,b) is labeled horizontal if |t(; — 6| = ei and is labeled 
vertical if — 6| =62- A dimer is labeled north (respectively south) if it is a horizontal edge 
with a black leftmost (resp. rightmost) vertex. A dimer is labeled east (respectively west) if 
it is vertical edge with a black leftmost (resp. rightmost) vertex. 

We are free to choose any Kasteleyn orientation and so we choose the same Kasteleyn 
orientation used in [26] which is given by the following: for notational purposes, let 

(2.9) Wi = {(xi,X2) G W : (xi +X2) mod 4 = 2i + 1}, 
and 

(2.10) B, = {(2/1,1/2) GW: (yi+y2) mod 4 = 2i + 1}. 

for i G {0, 1}. The edges incident to each vertex in Wq have edge weights —ia, —1, ia, 1 with 
the labeling given by proceeding clockwise from the edge leaving the white vertex in the 
direction of 62- The edges incident to each vertex in Wi have edge weights ia, 1, —ia, —1 with 
the labeling given by proceeding clockwise from the edge leaving the white vertex in the 
direction of 62 . These weights are depicted in Figure |3j 

Let An denote the vertices given in the above co-ordinates with side length '2n. Let diJ\.n 
denote the vertices in An H B of the form (0, 2k — 1) and drAn denote the vertices in An H B 
of the form (2n, 2k — 1) for 1 < k < n. 
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Figure 3. The complex weights with Wi e Wj and bj G Bj. 

For the rest of the paper, let K denote the Kasteleyn matrix for An with the above weights, 
that is, iiT : B X W ^ C where K(b, w) = K^^w for 6 G A„ n B and w G A„ n W with 

r (-1)'+^' if 6GBj,u; = 6+(-l)^ei,6^5^^„,Z G {0,1} 

(2.11) K{b,w) = l (-ly+J'ai iibeBj,w = b-{-l)^e2,b^diAn,le{Q,l} 

y otherwise. 

We can define the Kasteleyn operator iiT : — )■ for w G B and for / : B ^ M 

(2.12) Kf{v) = Y,K{w,u)f{v) 

where w ^ v means nearest neighbors. As An is finite, we shall interchange between operators 
and matrices without further mention. 

2.3. Height Function. The height function is defined on the faces of the square grid and 
is defined as follows: the height change between two adjacent faces is ±3 if there is a dimer 
covering the shared edge between the two faces and =Fl otherwise. As we traverse between 
two adjacent faces, we choose the sign convention to be +3(or respectively —3 ) if the left 
vertex is white (or respectively) black. The definition is consistent around each vertex, that 
is, the total height change around each vertex is zero. The dimer model is in bijection with 
the difference of heights between the faces of the graph. 

We denote the height function of the Aztec diamond by /i„ which is given by the following: 
the height change along the diagonal faces adjacent to the boundary of the Aztec diamond is 
±2. Suppose that the height function at (0, 0) is 0, then the height function at both (2r, 0) 
and (0, 2r) is equal to — 2r for 1 < r < n — 1 while the height function at both (2n, 2r) and 



8 



SUNIL CHHITA, KURT JOHANSSON, AND BENJAMIN YOUNG 



(2r, 2n) is equal to 2r — 2n for < r < n. For a face / in the Aztec diamond, we have 



m— 1 



(2.13) hn{f) = ^{-ir^H^x^m-i) 

i=0 

where 7 is a fixed path from (0,0) to / of length m, ^[i] is the i^^ index of the path and 
is the indicator function that a dimer is traversed from to 7[i + 1]. Figure [T] shows 
a realization of the height function for a random tiling of the Aztec diamond. 

2.4. Asymptotic Coordinates, average height function and the arctic ellipse. In 

this subsection, we introduce the coordinates of frozen region as well as the average height 
function and the southern domino process. As n — )• c«, for the (re-scaled) Aztec diamond 
with corners (0, 0), (1, 0), (0, 1) and (1,1), the boundary between the frozen and unfrozen 
regions is an ellipse [8] whose equation is given by 

(2.14) - + (Jl±l^ = 1 
^ ^ l-p p 

where u G [0, 1] is the horizontal coordinate, v G [0, 1] is the vertical coordinate, and p = 
1/(1 + a^). We let P C be the area bounded by the ellipse given in (|2.14l). 



The average height function was computed in [8j for a = 1 and can be extended for all 
values of a. From [8j, define 

11 V - (1-v) 

(2.15) P{x,y,p) = - + - arctan ^ ^ ^' 



2 vr — p2 _ — p)x'^ — py'^ 

which is the asymptotic placement probability of a north domino on an Aztec diamond 
with corners (1, 0), (0, 1), (— 1, 0) and (0,-1). The asymptotic placement probabilities of 
south, east and west dominos are given by P{—x, —y,p), P{—y,x, I — p) and P{y,x, I — p) 
respectively. Let 

(2.16) hr{x,y) = 2{yP{x,y,p)-yP{-x, -y,p)+{l-x)P{-y,x,l-p)+{l+x)P{y, -x, l-p). 

The asymptotic height function for an Aztec diamond with co-ordinates (0, 0), (1, 0), (0, 1) 
and (1,1) is given by 

(2.17) ~h{u,v) =~hri{l + u + v)/2,{l + u-v)/2) 

which can be checked by showing it verifies the two differential equations obtained from the 
expected heights across a face and its boundary conditions. 

We remark that J^J^/i(n, v) is equal to the so-called Archimedes measure [2j, which is the 
conjectured measure of the asymptotic density of I's of the permutation matrix of a random 
sorting network. More details on random sorting networks can be found in [2j. 

The boundary of the arctic ellipse is given by 

(2.18) V = 1 - u ±2y^{l - p)p{l - u)u + p{2u - I). 

The arctic ellipse can be parametrized by {u{k),v{k)), where v{k) = 1 — k'^u{k) and 

(2.19) u(k) = , 

^ ^ ^ ^ {l + a^){l + k^) + 2aVT+^k 

We will be interested in two parts of the boundary. The part where A; > will be called the 
nort/iern 6oMn(iarj/ and the part where A; E (—a~^{l+a'^)^^'^ , —a{l+a'^)~^^'^) = ( — 1/-^1 — p,—^/l — p) 
the southern boundary. 
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Figure 4. This figure shows the intersection of the limiting elUpse with 
the red dotted hne given hy v = 1 — k'^u. The south dominos on the blue 
dashed line represent the southern domino process. The northern boundary 
lies between A and B and the southern boundary between C and D 



See Figure [4] for an example of the northern and southern boundary and an explanation 
of the geometric meaning of k. 

2.5. Local Gibbs Measure. In this subsection, we briefly describe the local Gibbs mea- 
sures for the dimer model but a more formal account can be found in [29J and [28 ^ For all 
doubly periodic bipartite weighted dimer models embedded on the plane, [29j found that the 
dimer model is a Gibbs measure and also gave an explicit method to compute the entries of 
the inverse Kasteleyn matrix embedded in the plane. The results from [29j rely on using the 
smallest non-repeating of the graph unit called the fundamental domain. For the graph con- 
sidered in this paper, the fundamental domain has one black vertex and one white vertex. In 
order to describe the Gibbs measure, [29] introduced magnetic co-ordinates {Bx,By) where 
one increases the energy by e^"^ or e^^ if one crosses to the neighboring fundamental domain 
to the left or above. Conversely, if one crosses to the fundamental domain to the right or 
below, one decreases the energy by e^"" or e^^ . These magnetic co-ordinates are related to 
the average slope, that is, one can compute the Gibbs measures for different slopes, see [29]. 

We choose the fundamental domain of the graph in this paper (embedded on the plane) 
to be given by a white vertex, an edge parallel to 62 which has weight ai and its incident 
black vertex. The edges parallel to ei are given weight 1. This choice of weighting is the 



same as the choice of weighting in Section 2.2 up to Kasteleyn orientation. We denote 



the Gibbs measure of the model on this graph by ^a{Bx,By) where {Bx,By) is described 
above. Suppose that (2^i + 1, 2^2) is a white vertex and (2r/i, 2?72 + 1) is a black vertex for 
'ni^V2,(,i,C2 £ ^- Suppose that {Bx,By) are the magnetic co-ordinates as described above. 
Using techniques from one can find the entries of the inverse of the (infinite) Kasteleyn 
matrix, denoted by and they are given by 



1 f f r'7l-5l,,,'72-6 rjqnriz 

,2.20, A'-,(26 . 1. 26). (2,. 2,. pj,?^ 
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where P{z, w) is the so-called characteristic polynomial defined in [29j and is given by 
(2.21) P{z, w) = ai + + + aiw'^z^'^. 



One consequence of knowing (2.20) is that by using Theorem 2.2 we can compute the 
probability of any cylinder set. 

2.6. Gaussian Free Field. Another consequence of finding the inverse Kasteleyn matrix 
for domino tilings on the Aztec diamond and its corresponding asymptotics is that we are 
able to show that the fluctuations of the height function in the liquid region are governed 
by the Gaussian Free Field in the limit. This subsection introduces the Gaussian Free Field 
but a more in depth treatise can be found in [36J. 

The Gaussian Free Field on HI is a probability measure on the set generalized functions 
on M such that for any compactly supported test functions (pi, cj)2, {F, (pi) := F{z)(l)i{z)\dz\'^ 
is a real Gaussian random variable with mean zero and variance 



(2.22) K[{F,^i){F,cP2)]= [ Mzi)Mz2)G{zi,Z2)\dzi\^\dz 



2 



(2.23) G(zi,Z2) = -^log 



where 

' Zl- Z2 

zi-z^ 

For any finite sequence of compactly supported test functions {(pk}^ the Gaussian Free 
Field F oTiM. can be uniquely characterized [25J by its variance and that we can write 
(2.24) 

li " I Epairings ^F, <Pa(i)) {F, <^.(2)>] • • • </.,(„_i)> (F, <A,(^))] m is even, 

where the sum is over all (2m — 1)!! pairings 

2.7. Statement of Results. Let c : W x B — )• M denote the coupling function, that is 
C = K^^. The coupling function was originally computed in [16] for the uniform measure, 
by an explicit enumeration of a signed counting of domino tilings — the information encoded 
by the inverse Kasteleyn matrix. We generalize the coupling function so that we can consider 
different weights for vertical and horizontal tiles using the following method: By using the 
interlaced particle system [18^ IT9] , we guess the entries of the inverse Kasteleyn matrix and 
show that the guess is correct. We remark that the correlation from the interlaced particle 
system and the correlations kernel from the inverse Kasteleyn matrix are not in bijection 
for domino tilings of the Aztec diamond which differs to the lozenge tiling model |34j . We 
obtain 

Theorem 2.3. For x = (xi,X2) G W and y = (yi,y2) we have 



(2.25) cix,y) 
where 



h{x,y) forxi<yi + l 

fi{x,y) - f2ix,y) forxi>yi + l 



(2.26, M.y>-'-^^. II 

^yi/2 (^^^^yx2/2^^^_iY^n-X2)/2 
^ ' ^(-i+l)/2(y; _ z) (aw - l)(2n+l-|/2)/2(o + y;)(y2+l)/2^^ 
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and 

(2.28) f2ix,y)=^-^ „fe-^2-i)/2 / (i/a + z) 

^ ^ ^ ^' 2m (l/a + a + z)(2/2-^2+i)/2 

where £i is the positively oriented contour \z\ = e, £2 is the positively oriented contour 
\w — l/a\ = e and the contours do not intersect. 

The original approach for computing asymptotic quantities for domino tihngs on the Aztec 
diamond used a non-intersecting particle system [181 ES] adapted from the so-called DR 
lattice paths [32]. Although much insight and information is gathered from this approach, 
the method does not allow one to compute the local statistics of the underlying dominos 
exactly. 



For the proof of Theorem 2.3 we have a set of relations that the inverse of the Kasteleyn 
matrix should satisfy from the relation K.K~^ = I. Indeed a similar set of relations (without 
a boundary condition) was used in where the authors show that the correlation kernels 
of a family of interlaced particle systems satisfy the Gibbs property in the bulk. 



Using Theorems 2.2 and 2.3, we can now compute all local statistics of the dominos, in 
particular, we compute all local statistics of the southern domino process which is defined 
as follows: fix r, 1 < r < n. With a south domino on the line y = r we mean a domino with 
a white vertex w = (2x — 1, 2r) and a black vertex b = {2x, 2r + 1) for some x £ {1, . . . , n}, 
and we say that the south domino is located at x on the line y = r. The south dominos form 



a determinantal process by Theorem 2.2 and in particular so do the south dominos on the 
line y = r with the kernel given by the following lemma: 

Lemma 2.4. The kernel of the determinantal process given by the positions of the southbound 
dominos on a fixed line y = r in a random tiling of an Aztec diamond is 



(27ri)2 J^^ {aw - l)"-''(a + wy+^{w 

Proof. Theorem 2.3 gives the coupling function as ci{w, b) = fi{w, b) — /2(w, b), but in this 
case 

(2.30, /.K.)^<=^/ ''/;";'r' ..^o 

2TTt 1/a + a + z 

because the integrand is analytic at 2; = 0. The result now follows by applying Theo- 
The sign in front of the kernel will be — (— l)(^2~^i)/2 and we can remove the factor 



rem 



2.2 

{—lY^2-xi)/2 |-,y ^ conjugation. □ 

Consider a determinantal point process on a space A with correlation kernel K. Let 
< a < 1 and consider the point process obtained by removing each point in the process 
independently with probability 1 — a. We will call this new point process the thinned deter- 
minantal point process with correlation kernel K and parameter a. A way of modifying the 
original point process to obtain a point process with multiple points is by taking each point 
in the process 1^ independently with multiplicity m, where m is a geometric random variable 
with parameter (3, F[m = k] = {1 — /3)I3^~^, k > 1. We will call this (multiple) point process 
a thickened determinantal point process with correlation kernel K and parameter (3. 

Propostion 2.5. The thinned determinantal point process {xj} with correlation kernel K 
and parameter a is again a determinantal point process with correlation kernel aK, i.e. 

(2.31) E fe" = det(I - ct)aKlA), 
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for every ip G C^(A), = 1 — e ^ and A = supptp = suppcf). With the same notation the 
thickened determinantal point process with kernel K and parameter /3 is characterized by 



(2.32) 



E 



det / 



1-/3 



where {xj} is now the multi-set of points of the process. 

We will prove the proposition in Section [4} Note that the thickened process is no longer 
a determinantal point process since they are always simple point processes. 

We will now consider the southern domino process along the northern and southern bound- 
aries. Along the southern boundary the southern domino process is almost dense and we 
have to consider the dual process, the process of holes, instead. If we think of the locations 
in 1 , . . . , n of the southern dominos on the line y = r as positions of particles then the empty 
spaces, the holes, also form a determinantal point process with a kernel I — L, where L is 
the kernel for the particles. 

In our scaling limits we will obtain the Airy kernel point process which is a determinantal 
point process with kernel 

(2.33) 

Set 
(2.34) 



KAt{x, y) = / Ai(x + t)M{y + t)dt. 
Jo 



/3 



(2.35) 

and let A > be given by 
(2.36) - 



a 



1-/3 1 + a2 + akVlT 



"2 ' 



a{a + kVl + a^)^ 
(1 + a2)(aA;2 + A:\/rTa2)((l + a2)(l + P) + 2aA:\/r^^) 
with the plus sign for k > and the minus sign for A; < 0. 

Theorem 2.6. Let a > be fixed let A be given by (2l^ , a by (2l^ and /3 by l2l^ . 

Furthermore, let {xj} be the positions of the south dominos on the line y = {{1 — k'^u{k)ri\, 
where u{k) is given by (2.19). Consider a fixed k > 0. Then, the rescaled south domino 
process at the northern boundary, 

nu{k) 



AnV3 

converges to the thinned Airy kernel point process with parameter a. 



Next, consider a fixed k G 



(1 + a^Y'^, -a{l + a^y'-'^). Let {yj} be the positions of 



the holes in the south domino process, i.e. the dual south domino process, at the southern 
boundary. The rescaled point process 

yj — nu{k) 







Ani/3 



converges weakly to the thickened Airy kernel point process with parameter (3. 

Successive independent thinning and rescaling of a point process typically has a Poisson 
point process as its limit. If a tends to infinity, we see that the thinning parameter of the 
thinned Airy kernel a tends to zero. Hence we can expect that if we let a tend to infinity 
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with n (but not too fast) the south domino point process close to the northern boundary 
should converge to a Poisson process. This leads to the next theorem. 

Theorem 2.7. Fix A; > and let a = a{n), where a{n) — )• oo hut a{n) /ii}/^'^ — )• as n ^ oo. 
Set c{a) = 7r^/^(l + l/k)^/'^a'^/^ and let {xj} be the positions of the south dominos on the 
line y = [n(l — k'^uik))]. Then the rescaled point process 

nu{k) + c(a)nV3 _ ^. 

converges weakly to a Poisson process with density p{^) = \/ {1 — 

The condition on the allowed growth of a{n) is certainly not optimal but is an outcome 
of the proof. A similar result holds for the thinned Airy kernel point process on M: if the 
thinning parameter is sent to zero, the Airy kernel can be re-scaled to a Poisson point process 
on M which has a square root drop off. This result actually follows from the proof of Theorem 



2.7 Thus, we can think of the thinned Airy kernel point process as being an intermediate 
kernel between the Airy kernel and the Poisson point processes. We could also consider the 
behaviour of the leftmost south domino along the northern boundary. For a fixed a we will 
get convergence to the last particle distribution for the thinned Airy kernel point process, 
det(/ — Qi<'Airy)L2{^,oo)- When a goes to infinity with n but not too fast we will instead 
get one of the classical extreme value distributions in the limit, namely the last particle 
distribution in a Poisson process with density \/ {1 — ^)+. We will not give the technical 
details that are required to prove these results from the results above. 

For tiling models with frozen and unfrozen regions, the tiling measure is a local Gibbs 
measure. That is, by rescaling in a region so that the lattice remains but the boundary 
conditions are forgotten, the measure governing the dimer configurations is a Gibbs measure. 
For domino tilings of the Aztec diamond, we show that we can recover the local Gibbs 
measure in the interior of the unfrozen region. 

Theorem 2.8. Choose the rescaling so that the white vertices are given by 

(xi,X2) = {[ain] + 2^1 + l,[a2n] + 2^2) 
and the black vertices are given by 

(2/1,2/2) = ([ai"-] + 2r?i, [a2n] + 2r?2 + 1) 

forai,a2 G int(P) and for (^i, ^2, Vi^ ''12 £ ^- Then, the measure on domino tilings c onve rges 
to weakly to /ia(log \/ ai/(l — oi), log \/a2/{l — 02)) where is defined in Section 



2.5 



Our last result concerns the asymptotics of the fluctuations of the height function in the 
bulk. Following the well-established techniques for tiling models, we find that 

Theorem 2.9. For < a < 00 fixed, the fluctuations of the height function in the interior 
T> are given by the pullback of the Gaussian Free Field on H. 

The proof uses the fact that the moments of the height function can be found using the 
knowledge of the inverse Kastelyn matrix as used in [U [lOl ESI EU [27] . As a consequence, 
we have to compute the global asymptotics of the inverse Kasteleyn matrix. Tiling models 
with frozen and unfrozen regions, |3j gave an approach for proving the Gaussian free field 
for lozenge tilings models provided they are completely described by an interlaced particle 
system with densely packed initial configuration. Building on results from [34], [35] was able 



to prove universality of [3j. Much of the proof of Theorem 2.9 follows the approach taken in 
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[5] and so we omit some of the technical estimates as these can be found in [3l [35] . Other 
approaches for proving fluctuations of the height function arising from a tiling model have 
been considered in p[2j , where the author exploits a quadruple integral to bypass some rather 
technical estimates, and where the author considers the characteristic function of the 
height function using the Cauchy Riemann operators. The proofs of the results of 126\ |25] do 
not apply to domino tilings on the Aztec diamond due to the domino tilings studied there 
had the so-called Temperley boundary conditions. 



2.8. Overview^ of the Paper. The rest of the paper is organized as follows. In Section [3j 
we prove that Theorem 2.3 In Section |4j we give the the proofs of Theorem 2.6, Theorem 2.7 
and Theorem 2.8 The proof of Theorem 2.9 is given in Section [5j In Section [6j we give 
the relationship between the dimer correlation kernel and the particle correlation kernel, 
introduced in IT9l. 



3. Discrete Setting 



3.1. Some Comments on Theorem 2.3 Theorem 2.3 was first obtained in [16] for a = 1. 
The approach used was an explicit enumeration of 5)2"^""^^^/^ which counts the 
number of signed perfect matchings with w and b removed from the Aztec diamond where 
the sign comes from adapting the Kasteleyn orientation for the removed vertices so that the 
number of counter clockwise edges per face is odd (alternatively, the product of the edge 
weights around each face is negative). 

We can also give a version of the coupling function when the dimer model has weights 
Ti,r2,r2, and r4 for the north, east, south and west dimers. However, the partition function 
of the dimer model is function of the ratio of the weights ri/r^ and r2/r4 which means the 
extra weights are superfluous. 

3.2. Proof of Theorem ETH 

Proof of Theorem 2.3. The matrix is uniquely determined by the specific choice of the 



Kasteleyn matrix, K. We only need to check that KC = I where Cx,y = c{x,y). For 
6, y G B n An, we have 



(3.1) KC{b,y) = Y,Kib,x)Cix,y) 

= Y,K{b,x)C{x,y) 



where K{b,w) is given in (2.11). We have that the rows and columns of KC are indexed 
by black vertices. There are three cases to consider which are dependent on 6 = (61,62); 
< 61 < 2n, 61 = and 61 = 2n. In each of the three cases, we set j = (61 + 62 — l)/2. 



For < 61 < 2n, we have that each entry of ( |3.1| ) is given by 
(3.2) (-1)-'' (c(6 + ei,y) - c{b - ei,y) - aic{b + 62, y) + aic{b - 62, y)) 

For yi 7^ 61, we have that 



(3.3) (-I)-''' (c(6 + ei, y) - c(6 - ei, y) - aic{b + 62, y) + aic{b - 62, y)) = 
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which can be seen by directly manipulating the integrands of the functions fi{x,y) and 
f2{x,y), defined in ( |2.26D and ( |2.28D . For yi = 61, ([3^ is equal to 



(3.4) 
(3.5) 



bi+3b2 + 'iyi+y2 

2m 

_(_l)(j/i+?/2) b2 = y2 
otherwise 



2 (1 + a^)z 



£1 



62-2 / I \ 2 

— \- a + z ] az 
a 



by using Lemma 3.1 As (yi + 1/2) niod 4 is either 1 or 3, we have 

(3-6) KC{b,y) = ll n^l'^T'^^'Zl' ■ 

^ ^ ^ \ < 61 < 2n,2/i / 62 



For bi = 0, we have that each entry of ( |3.1[ ) is given by 
(3.7) (-1)^' (c(6 + ei, y)) + aic(6 - 62, y)) 



For 2 < yi < 2n, (3.7) is equal to 



(3.8) 



(-1) 4 



(27ri)2 

because the function is analytic at z = 0. 
For 61 = and yi = 0, (3.7) is equal to 

(3.9) 



(I + a ) (a + z)^2~ [az — I) 2 

. \ z ; ^—^dzdw = ^ 

82 J£i ^ ^ [aw — 1) 2 (a + 2 



3''2+!<2 
-1) 4 



(27r«)^ 



("2-1 2n — 1-62 

£2 -'fi ""^ ~ ^ (aio — 1) 2"^ (a + It;) 



(3.10) 



, ^ , 362 + y2 
(-1) ^ 

27ri 



2/2 -^2 -2 



1 + a^ z^T 



y2-'>2-2 / 1 



- + a + z 
a 



fc2-f2-2 
2 



The first integral in equation (3.9) is zero because the function is analytic at 2; = 0. By 
Lemma 3.1, the second integral in (3.9) is equal to — (— 1)^^ = 1 when 62 = y2 and zero 
otherwise. Therefore, we have 

(3.11) KC{h,y) 

For hi = 2n, we have that we have that each entry of (3.1 ) is given by 
(3.12) 



1 61 = 0, yi = 61 

61 = 0, 2 < yi < 2n 



(-1)^ (c(6 - ei, y)) - aic{b + 62, y)) . 



For bi = 2n and < yi < 2n — 2, (3.7) is equal to 
(3.13) 



6n+Sb2+yi+y2 

(-1) 



(27ri)2 



(1 + a^) w ^ [a + z) 2 (az — 1) 2 

I in ■y , ^ 27i+l-i;2 , , ^2 + 

e^JSi w-z ^^[aw-X) 2 [aj^yj) — 



-dz dw 



(3.14) 



+ 



6n+3b2+yi+y2 

(-1) 

2m 



Si 



tj2 —y2 ~ 

y2-b2-2 , r,^ K2-''2-2 / 1 \ 2 / 1 

a 2 (l + a'')^ 2 l--i-a + z\ \~ ~^ ^ 



2n — yi 
2 
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For the first integral in (3.13), we can compute the residue at z = using the following 
integral: 



(3.15) 



1 



{a + z) 2 (az — 1) 2 



-dz 



{a + w) 2 (^aw — 1) 2 



z^{w — z) w"- 

This follows from the fact that the sum of the exponents in the numerator is less than n and 
pushing the contour through infinity. Using (3.15) and the change of variables w = z + 1/ a, 
the first integral in (3.13) becomes 
(3.16) 



(-1) ^ 

27ri 



I'2~''2~2 

a 2 



(1 + a^)z 2 



+ a + z 



i'2-i/2-2 



+ z 



2n-yi 
2 



The above expression is exactly the negative of the second integral in (3.13). 



For bi = 2n and yi = 2n, (|3.12|) is equal to 

62-1 2n-l- 

(l + a^) w'^{a + z)^{az - 1) 2- 



(3.17) 



362 +K2 
(-1) 4 



(27ri) 



. 2ri + l — 2/2"' 

f2^£^l '^-Z z"'{aW-l) 2 (a + y;)~ 



Y^-z dw. 



Using (3.15), we can compute the residue at 2; = for ( |3.17| ) which gives 



(3.18) 



-1) ^ 
2m 



1 + a"')^ r — Trdw 

« ar« — 1 2 



1 &2 = y2 
otherwise 



using Lemma (3.1) with the change of variable w = 2; + 1/a. Therefore, we have 

(3.19) KC(h,y) = [ \ = o ■ 

^ ^ I 61 = 2n, < yi < 2n - 2 

Combining (Is^, (|3.1ll) and (|3.19|) shows that = C-i. 



Lemma 3.1. For /c G 



+ a + 2 



dz 



l+a2 





= 
otherwise 



Proof. For /c = 0, the left hand-side of (3.20) is equal to 

1 



(3.21) 



1 



dz 



2iTi J\^\=i z{l + a + z) l + a 



2 • 



dz. 



□ 



When k > 0, the integrand in (3.20) is analytic at z = and so the left-hand-side of (3.20) 
is zero. When A: < 0, we can move the contour to a small circle around — (a + 1/a) and use 
the fact that the integrand is analytic inside. □ 



4. ASYMPTOTICS OF DiMERS 



In this section, we will give the proofs of the results on the local asymptotics of the Aztec 
diamond. We start by proving Proposition 2.5 about thinned and thickened determinantal 
point processes. 
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Proof of Proposition \2. 5\ Let {yj} be the points of the determinantal point process with 



kernel K and let {nj} be independent Bernoulli random variables, 



1] = a. Let Ek 



denote the expectation for the determinantal point process and E„ the expectation with 
respect to the Bernoulli random variables. Consider the thinned process. Then, by Fubini's 
theorem. 



E 



Ea'E„ 



)) 



(4.1) 
(4.2) 



E 



K 



JJ(l-E„[l-e-"^'''(^^)]) 



E 



K 



det(/ - (jyaKlA) 



which proves (2.31). Next, let {mj} be independent geometric random variables, P[mj = 
k\ = (1 — (3)/3''~^, k > 1 and let Em denote the expectation with respect to these random 
variables. Then 



E 



E^Ek 



EKEra 



)) 



(4.3) 



since 



E 



K 



E 



(4.4) 



1 _ e-'^My^) = 1 _ (1 _ ^) ^ ^fc-ig-fc^{%) 

k=0 

1 _ g-'^fe) <P{yj) 



l_/3e-^fe) 1 -/3-/3</.(yj)' 



This proves (2.32). 



□ 



Proof of Theorem 2.6. By Lemma 2.4 the south domino process on the line y = r is a 



determinantal point process with kernel L given by (2.29). Let us first consider this process 
in a neighbourhood of the northern boundary, when r = [(1 — k'^u{k))n\, k > 0. (Below we 
will often neglect the integer part in this and in other expressions. It is not difficult to see 
that this is unimportant.) The kernel can be written 

Iff ^X2-nu{k) 

(4.5) L{xi,X2) = -j:—^ I '^^ I 



(27ri)^ ./ri Jv2 """" z^i-^C^) (a + w){w - z) 



^ng{z)-ng(w) 



where 

(4.6) g{z) = (1 - k^u^k)) log(a + z) + k^u{k) \og{az - 1) - u{k) log z. 

Here we choose the argument in the logarithms in the interval (0,27r). We see that when 



u{k) is given by (2.19), g{z) has a double zero at 

1 
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We can now use a saddle-point argument to analyze the relevant asymptotics of (4.5) and 
since this is a fairly standard Airy kernel asymptotics saddle point analysis we will not go 
into the details. For the integration contours in (4.5) we have chosen the steepest descent 
contours given by the level lines of the imaginary part of g{z) starting at Zc- It can be seen 
that we will have two ascending contours for the real part of g{z) which will leave in the 
directions e^'^*/^ and go to infinity. We will choose these for our i(;-integration. We will have 
two descending contours going from Zc to —a leaving in the directions e^^'^*/'^ and we use 
these for the z-integration. If we have the scaling 



xi = [nu{k) - 



X2 



(4.8) 
then 

(4.9) lim -Xn^/^z^^-'=^L{xi.X2) = a—^ [ dz [ 

(4.10) =aKM{^,v) 



[nu{k) - An^/^r/], 



dw- 



1 



iiz + w) 



uniformly for .^,77 in a compact subset of M, where a = Zc/{zc + a). Here T is given by 
z{t) = -te^'"-^)*, t < and z{t) = te*^, t > 0, with a fixed < 6* < 7r/3. Furthermore, A is 
given by 



(4.11) 



A = z,(-g(3)(^^)/2)i/3. 



A computation gives (2.36). To prove the result in Theorem 2.6 we observe that, with 
(j) = l-e-'l' and [n] = {1, . . . , n}. 



(4.12) E 



(4.13) 



■E,V'(«,) 



E 



n > 



nu{k) — Xj 
Ani/3 



m=0 



mi 



Xl,...,Xm&[n] j = l 



) Jmxm- 



Using the uniform convergence in (4.9) and Hadamard's inequality we see that (4.12) con- 
verges to 



(4 



" /'_1^m r ™ 

T. -^ / n<^(^^)d«t("^Ai(6.,0))„x™^"^ = det(/-./.ai^AilA) 



for every ^jJ E C^(M), where A = supp(/> = supp^. This proves the weak convergence 
claimed in the theorem. 

We turn now to the south domino process close to the southern boundary. Take r as 
before but with k G (— a~^(l + a^)^/^, — a(l + a^)~-^/^). We can write 



(4.15) 



L{xi,X2) 



1 



(27ri)^ 



dz 



7ri 



dw 



w""^ {a + zY{az-l) 



1 



7r2 



z^i {a + wY~^^ {aw — 1)^ ^ w — z^ 



where 'jr is \z\ = r positively oriented and < ri < r2, a < r2 < 1/a. Let a < < ri, 
r3 < 1/a and set 



(4.16) 



L{xi,X2) 



(27ri)2 



dz 



dw 



w 



(a + zY(az-l) 



1 



7r, 



z^i (a + wY^^{aw — 1)" ^ w — z 
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Here Zc is again given by (4.7). Note that now we have Zc G (—00,— a). If we move the 
ly-contour through the 2:-contour we obtain 



(4.17) z,^i-^^L(xi,X2) 



yXX-Xl 



(27ri) 



z 



X2—X1 



7ri 



z + a 



dz + L{xi,X2). 



The parameter (3 is given by /3 = —a/zc and we see that < /3 < 1. Set 

(4.18) M(xi,X2) = -/3"il(,,<,,), 

where \i^xt,<x2) is equal to one if xi < X2 and is equal to zero otherwise. Then (4.17) gives 

(4.19) .^1-^2, 

A correlation kernel L* for the dual point process is given by 
(4.20) 

If Xl,X2, 

(4.21) 



Z-^ -^L{xi,X2) = 6x1X2 - M(xi, X2) + L{xi,X2). 

point ; 

L*(xi,X2) = 5x1X2 - Z^'^''''^L{xi,X2) = M(xi,X2) - L(xi,X2) 

are the points in the dual point process we want to look at 

nu[k) " 



E 



n 



exp 



det(/ - f^g) 



where the Fredholm determinant is on £^({1, . . . , n}) and so is actually a determinant of 
a finite matrix. Here we have intr oduce d /(x) = 0((An^/'^)~^(x — nu{k))) and g{x) = 
lA((Ani/3)-^(x-nu(A;))). Now, by (4.20), we have 

(4.22) det(/ - fL*g) = det(/ - fL*) = det(/ - fM + fL) 



(4.23) 
(4.24) 



det(/ - fM) det(/ + (I - fM)-^fL) = det(/ + ^{fUyfL) 

j=0 



= det / + /^(M/)^L5 
since M is idempotent, M^'^^ = 0. Set 



(4.25) 



R = Y^iMfyi. 



Then, by (4.22), we have 
(4.26) 



det(/ - fL*g) = det(/ - f{R + L)g). 

N tt 

-^Ai(e,r?) 



In order to prove the result in Theorem 2.6 it suffices to show that, with the scaling (4.8), 
(4.27) Ani/3(i? + L)(xi,X2) - 



1 



uniformly for ^, 77 in a compact subset of M as n — t- 00. We will prove this under the assump- 
tion that tjj is also continuously differentiable, which suffices to show the weak convergence 
of the point process. 

We will make use of the following fact which again is proved by a saddle point argument 
very similar to the one discussed above. The only difference is in the choice of contours. As 
integration contours we will again choose level lines of the imaginary part of g{z) starting at 
Zc- It can be seen that we will have two ascending contours for the real part of g{z) which 
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will leave in the directions e^'^*/^ and go to 0. We will choose these for our w-integration. 
We will have two descending contours going from Zc to 1/a leaving in the directions e^^'^*/^ 
and we use these for t he z -integration. 



Let xi,rc2 be as in (|4^, r = [(1 - k^u(k))n\, k £ (-a(l + a^)"^/^, -a"i(l + a^)^/^) and 
a < < 1/a, < ri. Then, we have 
(4.28) 



- — ^ — / dz dwF(z) r — — r- 

(27ri)2 J^^^ J^^^ ^ ' z-i (a + wY+^aw - 1) 



"•w-z 1-/3 



uniformly for ^, r] in compacts as n — )■ oo, for F{z) = 1 or F{z) = (3zc/ [z — /3(1 — f{xi))Zf^, 
where D = F{zc)- 

By expanding M and rearranging the sum, we have 

n 

Xi,X2) 

J=l 
n 

= E /3"^"'^i(xi<m) • • • ife-i<.,)/(yi) • • • f{y,)L{y„x2) 

n n—xi 

Let ej{xi, . . . , xt-i) be the j:th elementary symmetric polynomial in t — 1 variables and write 
/ = (/(xi + 1), . . . , f{xi + t — 1)). By interchanging the sums, we obtain 

n—xi t—1 

R{xi,X2) = - p'f{xi+t)L{xi+t,X2)Y,^-iy-^e,^i{f) 

t=i 1=0 
4.29 ^ 

^ ^ n-xi t-1 

= - E + + ^2) 11(1 - /(xi + j)). 

t=i j=i 

Set 

(4.30) T{xi,X2) = -Yl /5*/(a;i)(l - /(a;i)*-iL(xi + t, X2). 

t=i 

We then have the following 

Claim 1. For all ^,r] in a compact subset o/M there is a constant C such that if we have 



the scaling {4-S) then 

(4.31) \Xn^/\R-T)ixi,X2)\< 



n 



1/3 ■ 
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We will prove the claim below. Using (4.16) we see that 
(4.32) 



t=i 



fc 

(2^i)2 

1 



dz 



7ri 



dw 



7r3 



z^i (a + wy~^^{aw — 1)"^'' w — z 



f{xi) z^^--^ 



1 - (27ri)2 
(a + z)''(az- 1) 



dz 

n—r 



dw 



7r3 



1-^(1- f{xi))Zc/z 



z^i (a + i(;)^+^(at(; — 1)" w — z 



-/(^i)- 



dz 



(2^i)2 

It now follows from ( 4.28[ ) that 
(4.33) Xn^/^T{xi,X2) 



dw 



w 



{a + zY{az-iy 



z - /3(1 - f{xi))zc -2^1 (a + wy+^{aw - l)"-*" u; - z' 
/3</'(e) 



1 



l-/3(l-,/.(0)l-/3 



^Ai(C,??) 



uniformly as n — )• oo. If we combine this with (4.31) and the fact that, again by (4.28) with 
F{z) = 1, that 



(4.34) 



An^/^L(xi, X2) 



1 



i^Ai(e,r?) 



we obtain (4.27), which is what we wanted to prove. 

It remains to prove the claim. We will use the following two facts. If ^, 77 belongs to a 
comp act subset, there is a constant i? > and a constant C such that if we have the scaling 
(0 then: (i) f{xi + j) = if J i [-Bn^/^ ,Bn^/'^] and (ii) \\n^/^L{xi + t,X2)\ < C for all 
t G [— Sn-'^/'^, That (i) holds follows immediately from the fact that (j) has compact 

support, and (ii) follows from (4.28) in the case F{z) = 1 since we have uniform convergence. 

Now, since (f) is continuously differentiable and has compact support. 



t-i 



(4.35) 



\l{l-f{x,+t))-\[{l-f{x{)) 

i=i i=i 

t-i f-i 

^(i-/(xi)^-i(/(xi+£)-/(xi)) n (i-/(xi+j)) 

t-l 

<Y^\f{x,+i)-f{x,))\ 



< 
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for some constant C. Thus, 

\\n^'^{R-T){xi,X2)\ 



(4.36) 



< p'\f{xi + t)\\Xn^'^L{xi + t,X2) 



t=i 



t-i t-i 

n(i-/(xi+t))-n(i-/(^i)) 

i=i i=i 



- „l/3 ^ - „l/3' 



t=l 



n 



since < ^ < 1. 



□ 



Above we have been concerned with the behavior of the south dominos at the boundary 
of the frozen region. We can also consider the behavior of the south dominos as we enter the 
bulk but still stay close to the boundary at a macroscopic scale. We then get as a scaling 
limit the thinned sine kernel point process with the same parameter a. We will not go into 
the details. 



Next we will give the proof of Theorem 2.7 which is concerned with the case when a grows 
with n but not too fast. In the case when a instead goes to zero with n but not to fast 
we will have convergence to the Airy kernel point process at the northern boundary for the 
south domino process. Note that the thinning parameter q— )-lasa— t-O. Ifa = 7/n for 
some fixed 7 > then we get instead the discrete Bessel kernel in the limit as n — t- 00. Again 
we will not go into the details. 



Proof of Theorem 2.1 . Set 
(4.37) 



L*(xi,X2) = ''^L{xx,X2), 



with L as in ( |45| ) and Zc given by ( [47| ). Take </> G Cc(M), < < 1 and let 

nu(k) + c(a)n^/^ — x\ 



(4.38) 



fix) = 4> 



c{a)n}l^ j 
for X € [n] = {1, . . . , n}. Furthermore, define 

(4.39) Mn{^, T]) = -c{a)n^/^Ll{[nu{k) - c{a)n^/^{C - 1)], [nu{k) - c{a)n^/^{r] - 1)]). 
Set 

nu{k) + c{a)n^/^ — (j + 1) nu{k) + c{a)n}/^ — j 



(4.40) 
for j € Z, and 
(4.41) 



c{a)n}l^ 



MO = ^ 



nu{k) + c(a)n^/^ — j) 
c(a)n^/3 
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for ^ € Ij, j £ Z. Let A be a compact subset of the real line such that suppcpn C A for all n. 
Then, 



(4.42) 



(4.43) 



E 



11(1 - 0(0)) 



E 



i-iy 

ml 



E Ylfixj)detiL*{xi,Xj))r 



m=0 xi,...,Xm&[n]j=l 

(-1)' 



E 



m 



m=0 

det(I-0„M„lA)L2 



n0n(O)det(M„(e„e,)) mxm 



(4.44) 

Assume that we can show that 

(4.45) W^nMnlAh^Q 

as n — )• oo, where || • ||2 is the Hilbert-Schmidt norm, and 

(4.46) M„(e,0^ 



uniformly for ^ € ^ as n — )• oo. The result then follows from (4.42) in the following way. 
if the operators Bn on L^(M) are trace class then the det2(/ + Bn) determinant and the 
Fredholm determinant det(/ + Bn) are related by 

(4.47) det(/ + Bn) = e*'^^"det2(/ + Bn). 

Also, if we have H-Bnlb — )• as n — )• oo, then we have det2(/ + Bn) — )• 1 as n — )• oo, see [TS]. 
Hence 

(4.48) det(/ - 4>nMnlA) = e*=^^"^"i^det2(/ - 0nM„U). 

It follows from (4.45) that det2(/ — (pnMn^A) — >• as n — t- oo, and from (4.46) that 

(4.49) tr0„M„U= / Mi)Mn{i,^)di^ [ mVi^-O+d^ 



as n oo. Thus, by ( |4.42[ ) and ( |4.48| ), 



(4.50) 



hm E 

n— >oo 



11(1 -'^(^^•)) 



which is what we wanted to prove. 

We turn now to the asymptotic analysis of M„ and the proof of (4.45) and (4.46). We 
will denote by C a generic constant that can depend on k and d in Claim [2] but not on n or 
a. Let A be given by ( 2.36| ) and define 

(4.51) M^^)(e,r/) = -Xn^^^L*{[nu{k) - An^/^^], [nu{k) - An^/^]) 

Claim 2. Let a be given by {2.3^ and fixd>0. Then 

(4.52) 

for aU^,r]e [-da^/^,da^/^]. 

We will prove the claim below. Set c(a) = c(a) / A and note that c 

(a) ~7r2/3a4/3(l + fc)2/3 

as a — )• oo. Then, we find that 

(4.53) M„(e + 1,7? + 1) = c(a)M^')(c(a)^,c(a)7?). 



\M^j^\C,r])-aKMi^,v)\ < ^ 
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Thus, with an appropriate fixed di > 0, we have 
(4.54) 

\\4>nMnlA\\2 = [ M^f Mn{x,yflA{y) dxdy < [ Mn{i + l,r] + lfdidr] 

JM? J[-(il,dl]2 



c a 



M^^\c{a)i,c{a)iif didii 



[-da4/3,rfa4/3]2 



where d = di-n'^/^{l + kf/^. Using (4.52) we see that 



(4.55) 



[-da4/3,da4/3]2 



Kki (e, V? d^dr] + 2(i2(a^/3a-2)2. 



Now a — )■ 0, since a = a{n) — )■ oo as n — t- c«, and we see that in order to prove (4.45) it 
remains to control the integral in (4.55). If we use the identity 

(4.56) 

and 
(4.57) 

we see that 



Kai {x, yf dy = Kai {x, x), 
Km {x, x) = Ai'ix)"^ - xAi {xf 



(4.58) 



< 



KAiitOd^ 

/oo 
Ai'(e)2-xAi iO'^dC 



[-da4/3,da4/3]2 



2a' 



2{da^'^fA\\-da^l^) + 2da^''' k\' {-da^'^f - Ai {-da^l^)M' {-da^'^) 



Now, as r — ?■ oo we have 
(4.59) 



Ai(-r) 



1 



„-l/4 



TT 



sin I -r^/^ + - I + . . . 



Ai'(-r) 



vr 



r^/'^cos f + ^ I + . . . 



and hence, since a ~ (1 + A;) a for large a, we obtain the bound 



(4.60) 



2a' 



-da4/3,da4/3]2 



Kai (C, V? d^dT] < Ca-^(da4/3)3/2 < Ca 



and since a = a{n) — )• oo we have proved (4.45). 
It follows from ( |4.52| ) and ( |4.53D that 



(4.61) 



|M„(C,e) - ac{a)KAi{c{a){i - l),c(a)(e - 1))| < Ca^'^-^ = Ca'^'^ 



for all ^ in a compact interval, and since a — )• oo as n — )• oo, (4.46) follows by using (4.57) 
and standard asymptotic formulas for the Airy function and its derivative. 

It remains to prove Claim [2] Let 7 > be given by ^^g^'^\zc) = —2 and recall that 
A = Zc/t- Let g{z) be defined by (|4.6[) and write 



(4.62) 



/€(C) = AnV3^1og [1 + 
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and 
(4.63) 



m) = MO + n{g{zc + iCn-^'^) - g{z,)). 



If we use (4.5) and make the change of variables z = Zc + '^Qn 
obtain 

(4.64) 



w 



Zc + "fion 



-1/3 



we 



a + Zc 



V 



{a + Zc + juin ^/^){uj - C) 



where C and P are the images of the steepest descent contours in theorem 
we have the estimate 



2.6 



If 



(4.65) 



n 2n 



n 



1/3 



n/2 



for ah sufficiently large n. C will lie inside \C,\ < a^n-*^/^, and using (4.65) and the estimates 
we describe below for the (^-integration, we see that we can replace P by the part of T) that 
lies inside |ci;| < a^n^^^. We denote this part by V also for simplicity. 

Let CI be the part of C in the disk |C| < n^/^^, the part in the annulus n^/^^ < ICI < n^^^^ 
and C3 the part in |C| > n^/^^. Let Cj and Cj be the parts of C* that lie in the upper and 
lower half plane respectively. We make the analogous definitions for V. We will consider 
estimates of F^{C,) on Cj, i = 1, 2, 3. The estimates on Cj are the same by symmetry and the 
estimates on V are analogous. The estimates that we need are 

(4.66) Jien{g{zc + lCn-^'^)-g{zc)) < -\\C? 
for aU C G Ci + C2 and 

(4.67) Ren(5(zc + iQu'^h - 5(2,)) <--n 

6 



7/15 



for all C £ C3 and n sufficiently large. To see this note that we can write 



(4.68) 
where 
(4.69) 



1 



n 



hliC) 



6nV3 



n 



1/3' 



Now, if we have |C| < n^^^^, then 
(4.70) 



|/ii(C)l < Ca^/^n^^/^lCr 



The estimate (4.70) follows from the fact that 

7s 



(4.71) 



n 



1/3- 



for \s\ < riJ since ^ <Ca If we have G Ci + C2, then = — ^Im i^^ + Im/ii(C) ^i^d if 
we write C, = re*^ this gives sin 30 = 3Im/ii(re*^) and hence, by (4.70), 

(4.72) I sin 36'| < Ca^l^rT^^^r < Cn-^^^. 
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Since Ci leaves Zc in the direction e^'^*/^ we must have cos 3^ > 2/3 for all large n. Thus 

1, 



(4.73) 
(4.74) 



Ren{g{zc + iQn ^1^) - g{zc)) 



-|Crcos30 + Re/ii(C) 
o 

<-^|C|^(l-Ca2/3n-V3|c|)<-i|C|3 



if C £ Ci + C2 and n is large. This proves (4.66). Since Yien{g{zc + ^C,n~'^^^) — g{zc)) is 
decreasing as we move along C in the upper hal f plan e starting at Zc we see that the value on 
C3 must be < —n^/'^^/Q by using the estimate (4.66) at the point where C meets |C| = n^^^^ 
(the endpoint of C2). This proves (4.67). 
We can write 



(4.75) 
where 

(4.76) 



h2{C) 



C-s 



1 



Xn^/^ Jo + 
and we see that if |C| < n'''^'^^, then 

(4.77) |/i2(C)l < Ca2/3n-V3|^|2_ 

We start by estimating F^{() on C3. If C ^ C3 then |C| < a^n^^^ and hence 



(4.78) 



C 



Ani/3 



< C7n^/^^ log 



n 



and combining this with (4.67) we see that Re-F^(C) < — j^n"^/^^ for large n and hence the 
contribution from C3 is negligible. 

Since \^\ < Ca^l^ and a < Cn^io we see that if C e C2, then Re/g(C) < Cn^/^^\C\. 
Combining this with (4.66) we see that ReF^(C) < — |C|^/12 for n large and hence the 
contribution from C2 is negligible. 

Thus, with an error that is much smaller than Ca~'^, we can replace Mn\£^, t]) by 



(4.79) M«(e,r/) 



a 



(27ri)2 



dC 



Ci+Ci 



dbJ 



a + Zc 



{a + Zc + jujn ^/^){cl! — C) 



Set S n = a {n) /n^/'^^ . By assumption 5„, — t- as n — )• 00. If |C| < n^/'^^ it follows from (4.70) 
and ( I^TtI that ^^(C) = -CV3 + + r„(C), where |r„(C)| < C5^r/^. Also 

a + Zc 



(4.80) 



{a + Zc + ^ojn 1/3) 



1 



if a; G X>i + X>i and thus we can approximate M„ (^, r/) with 



(4.81) 



(27ri)^ 



Ci+Ci 



a;-C 



Note that, if |C| > n^/^^, then ^|/|Cp < C a'^/^ /^2/i5 < (-.^4/3 ^^^^^ ^.^^ ^ negligible 
error, we can replace the expression in ( 4.81[ ) by a^Ai (C) '^)- This completes the proof of the 
claim and the theorem. □ 



We can now prove Theorem 2.8 
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Proof of Theorem 2.8. We show that the entries of the hiverse Kasteleyn matrix under the 



rescahng given in the theorem converge to the entries of the inverse Kasteleyn matrix given 



in (2.20). By using Theorem 2.2, we can conclude that the two measures agree on cylinder 
sets. 

Without loss of generality, we assume that xi < yi + 1 which means that K^^{x,y) = 
fi{x,y) for X = {xi,X2) and y = (yi,y2)- Let 

(4.82) gai,a2 (z) = 02 log(a + z) + {1 - 02) log{az - 1) - ai log z 

We have that g'^_^ (z) = for z equal to 

a^ai — ai + a2 + a^a2 — ± y^— 40^(1 — a\)a\ + (02 — ai + 0^(02 + ai — 1))^ 

^^•^^^ 2a(l - a2) • 

Suppose that Zc is the positive square root of the above equation. The square root term of 
Zc is given by 

— 4a^(l — a\)a\ + [0.2 — ai + 0^(02 + ai — 1))^ 
(4.84) /(ai-a2) 



\2 



a 



2 ^{a2^ax-\Y\a\\^a'')-a^ 



which is less than zero for ai, 02 G • This means that Zc is complex for ol\ and 02 in the 
interior unfrozen region. By a computation, we have that 



(4.85) |z,|=/-^=:ri. 
We also have 

' a + 



(4.86) 



azr — 1 



a2 

:j =: r-2. 

1-02 



We choose x = {xi,X2) and y = (yi,y2) as given in the theorem. As the double contour 
integrals are standard and the techniques of steepest descent are very well-established for 
this type of double contour integral formula, we omit the technical computational details for 
providing the bounds of the error terms. For fi(x,y), the contour £1 can be deformed to 
a contour from to Zc and a contour from Zc to which passes to the left of the origin. 
The contour £2 can be deformed to a contour which passes to the left of 1/a but also passes 
through the points z^ and Zc- This deformation picks up a single integral obtained from 
setting z = w in the expression from fi{x,y) while the remaining double integrals can be 
shown to be 0{n^^^'^) see [33] for example. We find that fi{x,y) is given by 

(4.87) : / z'l^-^^-'— ( ) z^^-^'^dz + 0(n-i/2 



27ri a + z \az — 1 



We can make the change of variables w = [a + z)/{az — 1). We find that fi{x,y) is now 
given by 

(4.88) / dw + 0(n-^/2 

where t{z) = (a + z)/[az — 1). The above integral converges uniformly to 

jY — l '|(5i+6+r)i+r?2)/2 [■ r- "6 ,,,'72-6 

(4.89) / / _l dzdw. 

{2my J\w\=r2 J\z\=ri a + w + z- awz 
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for ^1, ^2) ??2 £ ^- This is due to \a + w\/\aw - 
so one can use the residue formula to obtain (4. 

vc 

-6-1 ( 



1| < ri if and only w £ [t(zc),t{zc)] and 



variables z i— )• zai/i and w i— )• wa2/i in (4.89) and we obtain 



from (4.89). We can take the change of 



'2 



z w 



(2vri)2 7,^1=1 71.1=1 n^n,Pr2) 
Computing the probability of any cylinder set using the inverse Kasteleyn matrix given 



by (|2.20|) agrees with using the inverse Kasteleyn matrix given by (|4.90|) provided that 



ri and e^^ = r2. 



□ 



5. Gaussian Free Field 



In this section, we continue our study of the asymptotics of domino tilings of the Aztec 



diamond. We introduce the required complex structure for Theorem 2.9 We compute the 



necessary asymptotics of the inverse Kasteleyn entries and use them to prove Theorem 2.9 



The approach of the proof closely mirrors [3] and [35] but uses the determinantal process 
obtained from the inverse Kasteleyn matrix in the spirit of [27j. This is due to the height 
function of domino tilings of the Aztec diamond uses dominos as opposed to the interlaced 
particle system. To compute the asymptotics of the inverse Kasteleyn matrix in the bulk, 
we follow the approach given in [12] which computes the asymptotics of a similar double 
contour integral. As the proof showing that the height fluctuations is well-established by 
the methods from [3], we do not provide all the technical details of the proof - we refer 
the reader to [3] and [35] where the same type of estimates have been computed for lozenge 
tilings. 

5.1. Complex Structure. In this subsection, we introduce the complex structure as well 
as the notation that is required for the rest of the section. Define (7 : EI — )• C by 

(5.1) g{z\ (xi, X2)) = X2 log(a + z) + (1 - 3:2) log(az - 1) - xi log z. 

To every point (xi, X2), there exists at most one G IHI such that g'{Q\{xi, X2)) = 0. If this 
map exists, we can define 

(5.2) V = {{xi,X2)\^{xi,X2) £ M exists } 

We have the following the lemma for the map Q which follows from a similar calculation 
in [12]. 

Lemma 5.1. The map (xi, X2) 1— >• ^{xi, X2) is a diffeomorphism from D to H. Furthermore, 
we have T) = T). 



Proof. Since g'{^) = where g is defined in (5.1 ), we can split the equation into its real and 
imaginary components which give 



(5.3) 



a-\-\l 



Im 



1 




-Re 
-Im 



aQ-l 



The determinant of the above matrix is equal to alm(il)/|r2p(l/|a-|-ilp-|-l/|ar2 — Ip) which 
is non-zero because > 0. Therefore, we have 
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The second statement of the proposition is a direct consequence of the fact that the largest 
contribution from each integral when computing K~^{x, y) for x and y in the liquid region, 
comes from saddle point analysis which requires g'{z) = 0. 

□ 

For the rest of this section, we use the notation for / : 2? — t- M that 

(5.5) (</>,F) := / cl){x)F{x)dx. 

Jv 

where : P — )• K is a compactly supported test function. 

5.2. Asymptotic Entries. In this subsection, we compute the asymptotics of the inverse 
Kasteleyn matrix in the interior of the liquid region which is split into three lemmas com- 
puting the long range and short range asymptotic behavior for fi(x,y) and the asymptotic 
behavior for f2{x,y) for x = (xi,X2) and y = {yi,y2)- We also state without proof three 
lemmas that compute the asymptotics of fi{x,y) where one or possibly both x and y are 
near to or on the edge. The proof of these lemmas are technical arguments that can be found 
by mirroring the analogous lemmas in [3j and [35]. 

In order to state the main lemma of this subsection, we first find the correct deformation 
of the contours for the function g{z\{xi, X2))- The function g{z\(xi, X2)) has saddle points 
at the complex conjugates and z^. We can deform the contour £1 into a closed counter- 
clockwise oriented contour which passes through Zx and z^ and contains the origin. Label 
this contour Ti. We can deform the contour £2 into a contour which passes through Zx and 
and is asymptotically parallel to the positive real axis. Label this contour 
For the rest of the section, we drop the integer part from the computations. This has no 
effect on the computations. 

Lemma 5.2. For X = {Xi,X2),Y = (^1,12) with X,Y e mt{nV) and \X - Y\ > n^/^+^ 
where S > then 

(^•6) /^(^'^) = -^ E : — 7-r— + Oin-^^'^) 



27rn 



{0,1} (« + 4ni4^ - 4')^-9"{z%^)g"{Zy') 



where X/n x = {xi,X2) and Y/n y = (1/1,2/2) as n ^ 00, Zx solves g'{zx\x) = with 
lm{zx) > and 



(5.7) z: 



£x 



if£x = 
if£x = l 



and respectively for y. We choose \/ —g"{z^\x) so that the orientation of the line z = Zx + 
s / \/ — g" {z^\x) is tangent to Fi at Zx- We have a similar choice for ^ —g'\z}x\x). We choose 
\jg"{.Aj\y) -50 that the orientation of the line z = Zy + t/ ^g"{z^\y) is tangent to F2 at Zy. 

We have a similar choice for ^ —g"{Zy\y). 

Proof. The argument for proving the lemma uses the steepest descent method following a 
very similar approach to [3l Proposition 6.9] and \i2\ Lemma 6.1]. 
Recall that we are computing the asymptotic behavior of 

(5.8) fi{nx,ny) = / / — rdzdw 

{27ny Jg^ [w + a){w-z) 
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where nx = {nxi,nx2) and ny = {nyi,ny2). The contributions to fi{nx,ny) come from 
each of the four pairs of saddle points. By the above choice of contours and the choice of the 
orientations of g"{z^'^\x) and g'^Zy^ly), we have the contributions outside of the four saddle 
points are exponentially small. Depending on the choice of x and y, we either have that 



< 



The restriction |x — y| > n imposes that \zx 



> n 



-1/2+5 



For \zx\ < \zy\, the paths of steepest ascent and descent do not intersect. 

In this calculation, we only consider the contribution pair (zx, Zy) as the other three pairs 
follow from an analogous calculation. We introduce the change of variables 

t 



(5.9) 

and 
(5.10) 



Zx + 



w 



Zy + 



for t,s £ [ 

(5.11) ngiz) 
and 

(5.12) ng{w] 



n^/^, n"^/^]. We can expand the exponents in equation (5.8), namely, we have 

2 



f N , ng"{zx) 
ngizx) H ^ (z 



zxf + 0{z - zxf = ng{zx) -- + 0{n-^'^+^') 



ng{zy) + 



ng"{zy 



-{Z-Zyf+Oiz-Zyf 



ng(zy) + '- + 0{n-'/'^''). 



Zy\ > n 



-1/2+5 



by continuity, we have that \g" {zx\x)— g" {zy\y)\ > n 



-1/2+5 



Due to having | Zx 

We also have that w + a = Zy + a + 0{n~^/'^) and w 
uniformly on compact sets contained in intP, we have that the contribution from the saddle 
points Zx and Zy is given by 



Zy + 0{n~^/'^). Therefore, 



(5.13) 



1 



,4/2 



,4/2 



„n(g(2,|x)-g(2y|j/))+(-s2+t2)/2 



zdsdt + 0{n 



-l/2+35\ 



(27ri)2n 7„„4/2 y„„4/2 {zy + a){zx - Zy)yf -g" {zx\x)g" {zy\y)' 
By factoring out the constant terms, we can approximate the Gaussian integral which gives 

f,n{g{zx\x)-g{zy\y)) 



(5.14) 



+ 0{n 



-5/2N 



2TTn[zy + a){zx - Zy)y'-g"{zx\x)g"{zy\y) 

The other three contributions can be computed in exactly the same way. 

For \zx\ > \zy\, the contours Fi and F2 cross when applying the contour deformation to 



the paths of steepest ascent and descent. Due to the term l/{w — z) 
fi{x,y), in addition to the double integral we obtain an extra term 

I Qn{g{z\x)-g(z\y) 

27ri Jj z + a 



in the integrand of 



(5.15) 



-dz 



where ^ denotes the intersection of the deformation of the two contours Fi and F2 onto 
their paths of steepest descent and ascent. The double integral term can be computed 
asymptotically as given in the case \zx\ < \zy\ and the same expression holds. For the 
single integral given in (5.15), we wish to bound the integral and show that it is negligible. 
Manipulating g, we have that 

(5.16) Reigiz\x) - giz\y)) 



iVi - xi) log l^l - (y2 - X2) log 



az 



1 



a + z 



< 
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for < arg(z) < 7r/2 because \ {az — 1) / {a + z)\ < 1. We can choose a contour from ^ to ^ 
so that Rez < Re^, which gives 

(5,17) gnRe{g(z\x)-g{z\y)) ^ gnRe(a{5|x)-g(5|y)) ^ 

As lies on the contour of steepest ascent and descent, the right hand side decays exponen- 
tiahy with n. Therefore, we have 



(5.18) 



1 Qn{g{z\x)-g{z\y) 

-dz 



2-Ki Z + a 



< c- 



decays exponentiahy (compared to the asymptotics of the double integral). 

□ 

The next lemma gives estimates for /i {X, Y) when X and Y are close but still in the 
interior. 

Lemma 5.3. For X, y G int(nP) with \X — Y\ < n-^/^+^, we have 

^nKe(g{zx\x)-g{zy\y)) 



(5.19) \h{X,Y)\<C- 
where C is a constant and X/n — ?■ x and Y/n ^ y as n ^ oo. 

Proof. The proof of this lemma follows an argument from |121 Lemma 6.2]. Similar to the 
proof of Lemma 5.2 we have to deform the contours to pass through their appropriate saddle 
points and split into the cases when \zx\ > \zy\ and when \zx\ < \zy\. The exponential decay 



of the extra single integral from the case \zx\ > \zy\ holds from Lemma 5.2 It remains to 
bound the double integral, which is the same for \zx\ > \zy\ and \zx\ < \zy\, for each of the 
four pairs of saddle points. We only present the computation for the pair (z^, Zy). 

We apply the saddle point method from Lemma 5.2 using the same local co-ordinate 
transformation which gives the same Taylor expansions of g{z\x) and g{w\y). Contrary to 
Lemma 15.21 we have that 



s t 

(5.20) z-w = Zx- Zy+ = + = + O (n~^) . 

V-9"{zx\x)n y/-g"{zy\y)n 

where \zx — Zy\ < By continuity, we have that g"{zy\y) — g"{zx\x) is 0(n^^/^^"'^). 
Either \zx — Zy\ < rC^I"^ or \zx — Zy\ is of order 7i~i/2+a for < a < 5, which means that 
we either choose Zr — z,, or , ^ = H , * = when estimating z — w. In any case, 

y ^-g"{z^\x)n ^-g"{zy\y)n ^ ^ 

after bounding the integral we obtain a Gaussian integral which can be estimated to give 
the result. □ 

The next lemma finds bounds for f2{X, Y) which holds for all values of X and Y in the 
bulk and near the edge. 



Lemma 5.4. For X = {Xi,X2),Y = {Yi,Y2) with X,Y £ nV, we have 
(5.21) 1/2(^,1^)1 <C' 



gnRc(g(2c)) 



In 

where C is a constant, 

(5.22) g{z) = {y2 - X2) (log(az - 1) - log(a + z)) + {yi - xi) log z, 

Zc is chosen so that g'{zc) = with G H and X/n — )• x and Y/n ^ y as n —>• 00. 
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Proof. By performing the change of variables z i— )■ z — 1/a, we can rewrite the integral so 
that the contour of integration is around the origin. The upper bound follows by using a 
one-dimensional steepest descent argument around Zc and z^. 

□ 

The following lemmas consider the behavior of /i near or at the edge and finds bounds 
to limit the behavior. In order to state these lemmas, we distinguish the three regimes as in 
13 [35]: 

inixe&D \X — xn\ > n^/^ (X is in the liquid region) 

(5.23) cn^/^ < mfx^dv \X — xn\ < r?!"^ {X is close to the edge) 
inf^^g^X) l-'^ — xn\ < cri^/^ {X is at the edge) 

We do not include the proofs of these technical lemmas because they involve estimates of the 
double and single contour integrals /i and /2 which follow from very similar computations 
given in both [3, Proposition 6.11-6.13] and [35"1 Lemma 4.15-4.17]. 

Lemma 5.5. For X,Y ^ int(nP) with X/n — )• x and Y/n ^ y, with one or both x = (xi, 2:2) 
and y = (2/1, 2/2) in dV and assume that \zx — Zy\ is uniformly bounded away from zero in n. 
Then, we have 

„nRe(g{zx\x)-g{zy\x)) 

(5.24) \K-\X,Y)\ < C = 

n^/9"{zx\x)g"{zy\y) 

uniformly in n where C is a constant. 

Lemma 5.6. Suppose that X G nT> is close to the edge and Y € int(nD) and X/n — )• x = 
(xi,X2) G dV and Y ^ y = (^1,2/2) S int(P). Then, we have 

nKc{g(zx)-g{zy))+n^/'-^g{zx\X-x/n) 

(5.25) \K-\X,Y)\ < C 5 , , , , 

uniformly in n where C is a constant. 

It can be shown that n^/^g{zx\X — x/n) is order n^/'^ and is negative. 

Lemma 5.7. Suppose that X,Y ^ nT> are on the edge with X — )• x G dT> and Y ^ y £ dT> 
and that \zx — Zy\ is bounded away from uniformly in n. Then, we have 

nRc(g{z^)-g{zy))+n^/^g(z^\X-x/n)+n^/3g{zy\Y~y/n) 

(5.26) \K~\X,Y)\ < C 



2 
n3 



uniformly in n where C is a constant. 



5.3. Proof of Theoem 2.9, Before proving Theorem 2.9 we introduce some notation. Let 
T>n denote the unfrozen region of with lattice space 1/n. Let H^'^ be the set 

(5.27) {(xi, . . . ,Xm) G int(P™) : 3i,j such that |xj - Xj\ < n^^/^+'^j 

Proof of Theorem \2.9[ The aim is to show that we can write the moments of the pushforward 



of h under the map Q as a Gaussian with covariance given (2.23). For tiling models with 
frozen and unfrozen regions, this approach was first considered in [3] and later used in a 
universal setting in |35j . We follow the outline of these approaches and refer the reader to 
these papers for the full technical details. 

For X G T>n, let hn = hn{^~^{x)), where n is the size of An and /i„ is the height function 
of the Aztec diamond of size n. Let denote the average height of /i„ and let Hn = hn — h^. 
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Similarly, let /i^ denote the average height of hn, and Hn = hn — h'^- For compactly supported 
continuous test functions (pi G Cc(M^), we have 



(5.28) 



E 



i=l 



■ \dx 



2 

m I ) 



where the integral is the Riemann sum over the faces of that is, 

'Hr,..A = I Hnix)Hx)\dx\^ = Hn{x)(P{x). 



X face 
in D 



In order to estimate 1E[]^^-|^ Hn{xi)], we will first estimate lEiHi^i ^nixi)]. Using the defini- 
tion of the height function (2.13), we have the following: suppose 71, . . . ,7m are paths from 
distinct points on the boundary to xi, . . . where we assume that \xi — Xj\ > n~^/^"*"^, 
we have 



(5.29) E 



^H{x,) 



.i=l 



E [{X,, - E [X,,]) • • • (X,„ - E [X,J)] {-ly- 



H Hmom 



where ij £ jj denotes the ij^ index of the curve 7j. 



The above summand can be computed explicitly using the local statistics formula (2.6) 



and the corresponding formulas in ^26j. The summand (5.29) gives 



(5.30) 



i=i 

K-^{bi„Wi,) 



det 



K- 



where hi. and Wi- are the black and white vertices crossed by 7j. As we have chosen Xj| > 
^-1/2+5 ^ i.e. (xi, . . . ,x„) G \ T-C^^, we can write (5.30) using Lemmas 
highest order. Without loss of generality, assume that all 7j's only cross 



5.2 



to 



5.4 



up to 



lorizontal edges, 

and assume that the edges are of the form (Wq; Bq) or {Bi, Wi). Therefore, using the multi- 
linearity of the determinant and the above discussion on the entries, we have that the highest 
order term is equal to 



(5.31) 



E 



2m^_-^Yl+... + tm, 



det 







gn(g(42)-9(41)) 



gn(g(4»)-9{B/)) 



where z/ = zj. and Xj is the coordinate of bj, which is defined in Lemma 



5.2 
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One of the expansions of the above determinant gives 

om \ ) 



(5.32) 



^sgn(cr) 



•,en6{0,l} 



2^-nr=i5"(^?)n^: 



whereas the rest of the expansion gives a term, which we wih call G{z) which contains the 
oscillatory parts of (5.31). We can find G{z) by subtracting (5.31) by (5.32) 

Therefore, we have an expression of E[J3if(xj)] and moreover, we have an expression of 
E[J3i/(xi)]. Using the convergence of the Riemann sum and writing 

1 1 



dz- 



2,g"{^{z^))n 



for the map Zi i— t- ^{zi)^ we have 



E 



(5.33) 



^ sgn(f7) 



E 

o-eSn ei,...,e„e[0,l] 



(27r)^ 



+ 0(n-i). 



by using the Riemann Lebesgue theorem which gives the integrals containing G{z) converge 
to zero as n tends to infinity. 

Using [5l Lemma 7], which says that 



(5.34) 



E n 

cr£m— cycle i=l 



1 



0, 



we have that the only terms remaining from (5.33) are the fixed point free involutions. We 
obtain 



E 



m/2 

pairings a i=l 



where G will be found when we compute the variance. 



There are two technical considerations for E 



>5^(a;a,+i)) 



one where some of the points 



Xi are near to the edge and the second when some of the points Xi are close together, i.e. 
Xi E K^^. In either case, we outline the main points of the argument. For a full technical 
proof of these considerations, one can mirror the approach in [3] or |35j . 

For the case where some of the xis are close to the edge, one can approximate E HI^i H{xi) 
using the computation given above in the following way: as we are free to choose the paths 7^ 
and their boundary end points, we make the choice so that the boundary endpoints of each 
7i are closest to Xi if Xi is close or on the edge. For example, if nxi is n^/^ from the ed ge, then 
the path length of 7^ (before rescaling) is order n^^'^. Using the estimates in Lemma 



and the appropriate choice of each 7j , then the approximation of E 



5.5 



5.7 



using the 



Riemann sum as above will tend to zero (due to summing over too fewer points). 
When some of the Xj's are in V.^^, we can approximate the E YYiLi H{xi) 



by first 



considering all the Xj's not in Ti"!^^ using the approach given above and then add an approx- 
imation in the neighborhood of Xi for Xi in 'HJfr . The latter approximation can be bounded 
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above by using Lemma 5.3 and Lemma 5.4, Thus, one is able to have an upper bound of 
E YliLi H{xi) by where e > 0. More formal versions of this argument can be found in 



[31 [271 E5|- A consequence of this argument is that integrating E HI^i ^i^i) 
tends to zero as n tends to oo. That is, 



over Til 



(5.35) hm E / rr(/),(xi)i?„(xi)|dxi| 

Finahy, we have 



0. 



,S i=l 



hm E 

n—^oo 



.i=l 



hm / h / 



hm 



m rn/2 

V / T\(t>{xi)T\G{x^(^i),X^(^i+i))\dxi\^---\dXn\ 

— Jt>'^ 

lirings " j=l i=l 



pairing! 



where • • • represents the integrand given in equation (5.28). It remains to compute the 

1 



variance. Using Equation (5.33), we have that 

m 

\{H{x, 



E 



.1=1 



(72) [Zi - Z2Y 



1 



(2vr)- 



In 



r2(xi) — r2(x2) 



r2(xi) — ^(x2) 



2tt 



G{n{xi),n{x2)). 



This means that G = Gr as required. 



□ 



6. DiMERS TO PARTICLES: GUESSING K ^ 



It is easy to derive a correlation kernel for the zig-zag particle point process in the Aztec 
diamond from the knowledge of the inverse Kasteleyn kernel using Theorem |2.3[ The trans- 
formation between the coordinate system used in [TU] for the particles (ui,U2) and the 
coordinate system {xi,X2) used for the dimers in this paper is 



(6.1) 



Ul = X2 

U2 = {X2 - Xi + l)/2. 



The correlation kernel for the zig-zag particles is given by |19| 

(6.2) Kniui,U2;Vi,V2) = Kniui,U2;Vi,V2) - (l)m,v:,{u2,V2) 

where 
(6.3) 

1 r dw r dz z""^ (1 -az)"-'^+^2(l + a/z)* w 
w 



kn{2r - £i,U2;2s - £2,^2) 



z (1 _ au;)"-''+=i (1 + a/wY w - z 
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(6.4) 



')2r-ei,2s-e2{U2,V2) 



l(2r-gi < 2s-£2) 

27H 



71 



(1 - az)'-'^+^2-ei dz 
(1 +a/zy-'' 7 



and ei, £2 £ {0, 1}, a < ri < 1/a, < r2 < n. 

Note that we have a zig-zag particle at 6 G B if and only if a dimer covers the edge (5— ei, b) 
or the edge (5 — Cj, 6). Thus, for bi, 62 G B it follows from Theorem 2.3 that 

(6.5) 

2 

P[particles at 61 and ^2] = K{bi, bi — er^)K{b2, ^2 — ^r2) det(K^^(6j — 6^^, bj))ij=i^2 

ri,r2=l 



det I ^K(6j,6j - er)-fir ^(ftj- 6^,6^)1 

\''=l / i,j=l,2 



Let 61 have particle coordinates (ui,U2) and dimer coordinates {xi,X2) and let 62 have 
particle coordinates (^1,^2) and dimer coordinates (^1,^2)- We also have 

Kn{ui,U2;Ui,U2) Kn{vi,V2;Ui,U2) 
Kn{ui,U2;Vi,V2) Kn{vi,V2;Vi,V2) 



(6.6) P[particles at 61 and 62] = det 



(or the transpose of the same matrix) . 

It is not obvious how to make an identification between the elements in the determinants 
in (6.5) and (6.6), but it reasonable that something like 
(6.7) 

^^^^Y-^^^Kn{vi,V2; ui, U2) = K{bi,bi - ei)K~^{bi - ei, ^2) + K{bi,bi - e2)K"^(&i - 62, 62) 

1p[Ui,U2) 

could hold with an appropriate 'i/'- Inserting the coordinates and using (6.1) we get 

^^(2/1, 2/2) f y2 - yi + 1 X2-X1 + 1 

-J^n 2/2, ;X2, 



V'(xi,X2) 

(6.8) 

= _(_l)-(-i+-2-i)/2 (K-i((^, - 1,X2 - 1), (yi,y2)) - aiK-\{xi + 1,X2 - 1), (yi,y2))) 

Note that for the zig-zag particles on the black vertices, we have xi G 2Z and X2 G 2Z + 1. 
Hence, we can write ui = X2 = 2s — 1, vi = y2 = 2r — 1, i.e. r = [y2 + l)/2, s = (x2 + l)/2. 
It follows from (6.3) that 

/ 2/2-2/1 + 1 a;2 - xi + 1 
2/2, ;2;2, 



(6-9) = ^ ..2 

We now use the identity 



r ?i;2^i/2 (az- 1)"-(^2+i)/2+1(^ + q)(x2+i)/2 ^ 

7, '^'^y^, {aw - l)n-fe+l)/2+l(y; + a){2'2+l)/2 ^{w - z) 



(6.10) 



z + a 
W2 



(^_l)-(^i+^2-l)/2^_-|^j(a::i+a::2)/4-(yi + 



3/2 )/4 



(_l)(a^l-l+a^2-l+yi+i/2)/4' 
^(^i-l-l)/2 

( _ 1 ) {^1 +i+^'2 - +3/2 ) /4 

^(x-i+l-l)/2 , 



in (6.9). Note that we have 

[lo.ll) (^_l-^{y2-X2) /2 ^_Y^{xi+X2) / 'i-{yi+y2) / i _ (-_]^^(3::i-a;2)/4-(yi-j/2)/4_ 
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(6.12) 



'^{xi,X2) 



Kn y2, 



y2 - yi + 1 X2- xi + l 
-i: ; X2, 



(6.13) = (-i)-(-i+-2-i)/2 (f,(^(x, - 1,X2 - 1), (yi,y2)) - ai/i((xi + 1,X2 - 1), (yi,y2))) 



where /i is given by (2.26). If we similarly write 
V'(yi,2/2) 



tp{xi,X2) 



?/2 - yi + 1 a;2 - xi + 1 



(6.14) = (_i)-(-i+-2-i)/2 (/^((xi - 1,X2 - 1), {yi,y2)) - ai/2((x2 + 1,X2 - 1), (^1,^2))) 

using (6.4), ( |6.10| ) and (6.11) we get 

(6.15) 



('_l)(2;i+a^2+J/i+J/2)/4 

/2((2;i,x2), (yi,y2)) = 7^. l(y2 < ^2) 



27ri 



71 



(z + a)fe-^2+i)/2 z ■ 



That this is the same as the /2((xi,X2), (?/i, J/2)) defined by (2.28) follows from the following 
computation: 



1(2/2 < X2 



(z + a)fe-^2+i)/2 z 



(z + a)fe-^2+i)/2 z 
l(y2 < X2)l{yi < xi) / z 



(,.-.,+i)/2) (az-l)fa^-^^i)/^ 



71 



-l(y2 < a;2)l(yi < xi) 
-l(yi < xi 



'|2— l/a|=£ 

-l(yi < xi)a(^2-^2-i)/2 



(z + a)(2'2-^2+i)/2 2; 
(,,-x,+i)/2) (az-l)fa^-"^-i)/^ 

|.-l/a|=e (2 + a)(j/2-X2+l)/2 ^ 

(,,-..,+i)/2) (az-l)fe-"^-i)/^ 

(z + a)(j'2-^2+i)/2 2 

{z + l/a)(^i"^i"^)/2)z(^2-a;2-i)/2 
|^|=e (z + a+ l/a)(f2-^'2+i)/2 



dz. 



If we combine ( |6.2| ), (6.8), ( |6.13 ) and (6.14) we are led to the conjecture 
(6.16) Er"^((xi,x2), (yi,y2)) = /i((2;i,x2), (yi,y2)) - /2((a;i,a;2), (^1,^2)) 



which is shown to be true in the proof of Theorem |2.3[ 

There is another approach based on looking instead at the probability of seeing a particle 
at a specified white vertex and a particle at a specified black vertex. This leads to the 
following relation between the particle kernel and the inverse Kasteleyn matrix, 
(6.17) 

y2-yi + l X2-X1 + 1 

2 2 



K-i((xi,X2),(yi,y2)) = -(-l)(-i--2+yi-y2+2)/4^^(y^^ 



where Kn is given by (6.2). Compare the relation found between the particle kernel and the 
inverse Kasteleyn matrix in [34j. 
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